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We examine the variational and conformal structures of higher order theories of gravity which 
are derived from a metric-connection Lagrangian that is an arbitrary function of the curvature 
invariants. We show that the constrained first order formalism when applied to these theories 
may lead consistently to a new method of reduction of order of the associated field equations. We 
show that the similarity of the field equations which are derived from appropriate actions via this 
formalism to those produced by Hilbert varying purely metric Lagrangians is not merely formal 
but is implied by the diffeomorphism covariant property of the associated Lagrangians. We prove 
that the conformal equivalence theorem of these theories with general relativity plus a scalar field, 
holds in the extended framework of Weyl geometry with the same forms of field and self-interacting 
potential but, in addition, there is a new 'source term' which plays the role of a stress. We point 
out how these results may be further exploited and address a number of new issues that arise from 
this analysis. 
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I. INTRODUCTION 



According to the standard variational principle (Hilbert variation) that leads to the Einstein field equations of 
general relativity, the gravitational action J R \J—g is varied with respect to the metric tensor of a spacetime manifold 
that is taken to be a four-dimensional Lorentz manifold (A^,g) with metric g and the Levi-Civita connection V. 
However, in many instances (see Ref. |l| for a complete review) one considers a Lorentz manifold with an arbitrary 
connection that is incompatible with the metric, i.e. 7^ 0. A motivation for such a generalization was initially 
inspired by early work of Weyl (Ref. ^). In this case, one considers variation of an appropriate action with respect 
to both the metric components gab and the connection coefficients F^^ without imposing from the beginning that 
^hc usual Christoffel symbols. In the current literature this variational principle, where the metric and the 

connection are considered as independent variables, is referred to as the first order or metric- connection or simply 
Palatini variation. (For an historical commentary of this principle of variation and related issues, we refer to Ref. |.) 

Such alternative variational methods were first analysed in the framework of nonlinear gravitational Lagrangians 
by Weyl (Ref. |), Eddington (Ref. |) and others (Ref. |- |^). In an effort to obtain second order field equations 
different from Einstein's, Stephenson (Ref. ||,||) and Higgs (Ref. 0) applied the first order formalism to the quadratic 
Lagrangians R^, RabR"^'^, RatcdR"''"^'^ and Yang (Ref. |) investigated the Lagrangian RabcdR'^'"^'^, by analogy with 
the Yang-Mills Lagrangian. However, Buchdahl (Ref. |) pointed out a difficulty associated with this version of the 
metric-connection variation which is related to imposing the metricity condition, i.e. the connection coefficients equal 
to the Christoffel symbols, after completing the variation, and subsequently constructed specific examples showing 
that this version of the first order formalism is not a reliable method in general (Ref. |l^). Van den Bergh (Ref. [Tl| ) 
arrived at a similar conclusion in the context of general scalar-tensor theories. The R + aR^ theory including matter 
was investigated in this framework by Shahid-Saless (Ref. |l^) and generalized to the / (i?) case by Hamity and Barraco 
(Ref. ^3|). These authors also studied conservation laws and the weak field limit of the resulting equations. More 
recently, Ferraris et al. (Ref. |l^ showed that the first order formalism applied to general / {R) vacuum Lagrangians 
leads to a series of Einstein spaces with cosmological constants determined by the explicit form of the function /. 
Similar results were obtained in the case of / (i?ic^) theories by Borowiec et al. (Ref. [l5| ). 

A consistent way to consider independent variations of the metric and connection in the context of Riemannian 
geometry is to add a compatibility condition between the metric and the connection as a constraint with Lagrange 
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multipliers. In vacuum general relativity, this constrained first order formalism results in the Lagrange multipliers van- 
ishing identically as a consequence of the field equations (Ref . |l6| ) . This method was applied to quadratic Lagrangians 
with the aim of developing a Hamiltonian formulation for these theories in Ref. ^ 

Consider a Lorentzian manifold (A^,g, ^) of dimension D where ^ is an arbitrary symmetric connection. Hence, 
/Vg =/= that is, the connection coefficients are functions independent of the metric components and the Ricci tensor 
is a function of the connection only. In the case of general relativity without matter fields, varying the corresponding 
action 

S^jh^dPx, L = g™i?„„, (1) 

one arrives at the well-known result that variation with respect to the metric produces vacuum Einstein's equations 
whereas variation with respect to the connection reveals that the connection is necessarily the Levi-Civita connection 
(provided that D ^ 2). The integral (|l|) is taken over a compact region U of the spacetime {M, g, 7^) and we assume 
that the metric and the connection are held constant on the boundary oiU. In the sequel we omit the symbol (Px 
under the integral sign and set w g. Gothic characters denote tensor densities, for example Qab '■— wgat- 

In the presence of matter fields there is an ambiguity because the compatibility condition between the metric and 
the connection does not hold. The matter Lagrangian depends primarily on the field variables ip, and assumes a form 
that is a generalization of its special relativistic form, which is achieved via the strong principle of equivalence and 
the principle of minimal coupling, according to the scheme r]ab — *■ gab and (9 — > (the order of the two steps being 
irrelevant as long as the connection is the Levi-Civita one). Variation of the total action 



S = J [«(,g,r)+£,„(.g,^,^V')], (2) 

gives the following pair of equations 

Gab^Tab--= T-^' 



SlNdQ'^ + SlPaQ'" - '^NcQ^' - (3b) 

ab 

These equations are inconsistent in general unless the matter Lagrangian does not depend explicitly on the con- 
nection, i.e. (5£m/<5r^j = 0. (It is interesting to note that in the case D = 2 in vacuum, since the equation 
^Ib = {aJ + 1/2 {^IQb + 5lQa - QabQ"), with Qa := - p a'^^w ^ -^alnw + Fa, and Va := T^ab-. lias a vanishing trace, 
(1 — D/2) {da^T^\/—g — Fa) = 0, the Ta part of the connection is undetermined (Ref. |l8|).) 

In Section ||, after a brief review of the unconstrained metric-connection variational results in higher order gravity 
theories, we show how the constrained first order formalism is used to prove that the field equations of these theories 
can be given in a reduced form that makes the comparison to the usual Hilbert equations direct and point out that 



such a correspondence is due to the diffeomorphism covariance property of the associated Lagrangians. In Section III , 
we prove that in Weyl geometry, higher order gravity theories are conformally equivalent to general relativity plus a 
scalar field matter source with the same self-interacting potential as in the standard Riemannian case, together with 
a new 'source term' that arises as a result of the presence of the Weyl covariant vectorfield. In the last section, we 
comment on the usefulness of the extended form of the conformal equivalence and point out how these results may 
be further exploited to develop this framework. 

II. CONSTRAINED AND UNCONSTRAINED VARIATIONS 

We begin with a Lagrangian that is a smooth function of the scalar curvature R and vary the corresponding action 

S = I wf {R) (4) 



with respect to the metric tensor and the connection to obtain respectively 



f'R{ab) - ^fgab = 0, (5a) 
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Fa H'g"^) = 0. (5b) 



Explicitly, the F— equation (5b) reads 

(da lnw + (In /')' daR-Ta) gbc - dagbc + rZ9mc + TZgmb = 0, (6) 
and so we can solve for Fq and substitute back in (H) to find 

dagbc = ^Ta9mc + ^Ta9mb, (7) 

where we have introduced a new metric Tjab ■— f dab, with conformal factor /'. This means that F is the Levi-Civita 
connection for the metric g. 

Equation (|5a| ) is more straightforward. On the one hand, its trace /' {R) R — 2f (R) is satisfied identically if 
f (R) = ai?^~^(up to a constant rescaling factor a) and so (^ ) becomes Rab ~ (l/4)i?.gafc — 0, (provided that 
/' (i?) ^ 0) which finally gives Rab — (l/2a)5ah = so that the underlined manifold is an Einstein space with constant 
scalar curvature R = g°'^Rab = 2/a. On the other hand, one could regard the above trace as an algebraic equation 
in R with roots pi, p2, ■ ■ ■■ This situation was analysed by Ferraris et al. (Ref. ^4|) who showed that such an analysis 
leads to a series of Einstein spaces each having a constant scalar curvature (see also Ref. |lO|) . 

By a completely analogous procedure one finds Einstein spaces for the choice L = / (r) where r — QabQ"^'', and Qab 
is the symmetric part of the Ricci tensor and also for the Lagrangian L = f (K) where K = Rabcd.R"''"^'^ ■ Note that in 
this last case, varying the corresponding action with respect to the metric and the connection, one obtains 

~ 2j9ab — J i^a ribklm + J ^ alm^kb + ^/ lam^k 6 — t) (S) 

and 

(wfRj"^^') = 0, (9) 

with trace f (K) K = f (K). Hence, either f [K] — aK identically or, given a function /, the trace is solved 
algebraically for K . In contrast to the previous cases, there exists no natural way to derive a metric g from the field 
equation (^ unless the Weyl tensor vanishes (Ref. |l^). Notice that the field equations derived from the Lagrangian 
R[ab]R'^'' by the Palatini variation impose only four conditions upon the forty connection coefhcients and leave the 
metric components entirely undetermined (Ref. |lO[). 

Let us now introduce a vectorfield Qc called the Weyl covariant vectorfield and assume a linear metric-connection 
relation 

Fc9ab = ~Qc9ab, (10) 

and define 

Cat — ^ab ^ {ab} i (H) 

where, 

Kb - {lb} + ^5'" {Qb9am + Qa9mb " Q^9ab) ■ (12) 

The constrained first order formalism consists of adding to the original Lagrangian the following term as a constraint 
(with Lagrange multipliers A) 

Lc {g, r, A) = A^"*" [Knn ~ {mn} ~ ^mn] ■ (13) 

For instance, in Riemannian geometry ( |l3|) takes the form 

Le (5,r,A) = A- [r;,„-{LJ], 

while, in Weyl geometry 

1 



ic (ff,F,A) = A, 



^mn {mn} {Qngms + Qmgsn Qsgmn) 
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As an example, consider any of the previous test Lagrangians L (g, T, tp) and vary the resulting action 

J w[Lig,T,^) + L,ig,T,A)], (14) 

with respect to the independent fields g, T, A and ip. We find the 5— equations 

6{wL) 



Sg 



ah 



+ wBab = 0, (15) 

r 



where Bab is defined by 



and the F— equations 



Bab ■= — 2 ^™ [Abam + Aamb — Amah] , (16) 



5L 



ab 



+ Af = 0. (17) 
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Variation with respect to the matter fields ip yields their respective equations of motion. 

In the Riemannian case in particular we have (some of the expressions given below appear incorrectly in Rcf. p^ ), 
for the Lagrangian L = R^, 

^R^g"'' ~ 2RR''^ + B"'' = (18a) 

A,"^ = (25'^"^ - g'^-^^l ~ g^^Ht) JQ ,,,R (18b) 

B"'' = -2g'''' pR^lp^p^R, (18c) 

for the Lagrangian L — RmnR"'^", 

^i?„„i?™"(?'^^ - i?'^"i?''„ - RJR""" + B"'' = (19a) 

A,^' = 2]^,R^' - <^:^n^i?"' - 'Jc^mi?"" (19b) 

^ab ^ _pj^ab ^^^^^bj^am _ ^ab ^ ^ ^j^mn ^ ^^g^) 

for the Lagrangian L ~ RmnrsR™""^^^ ^ 

ii?™„,,i?"""-^g'"' - 2i?''""'-i?''„„,, + S-^" = (20a) 

A,''^ = 2 ^^i?,""" + 2 ^mi?,^"" (20b) 

fi*^*- = 4/7„;^„i?"'"''", (20c) 



and for the Lagrangian L = f (R), 



\fgab - fR{ab) +Bab=0 (21a) 



(25''-C - /-<S^ - g'^-^51) = A,''^ (21b) 
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gab ^ _gabpj^>^ 77'' /7°/'. 



(21c) 



It is straightforward to obtain the correspondence with the Hilbert case by substituting the BabS in the first equation 
in each of these above cases. They read respectively 

ii?V''--R^'''' + V''V"i?-5°^ni? = (22a) 
]^R^nR""'g°''' - 2i?''""'"i?^„ + V''V°i? - □ii'^'' - ]^nRg°-^ = (22b) 
ii?^„,,i?'""'-^<?'^'' - 2i?"'""i?„„; + 4V„V„i?'^""'' = (22c) 

f'R(ab) - \f9ab - VaVbf + gab^f = 0. (22d) 

Strictly speaking, compared to the usual Hilbert variation, in all cases considered so far using the first order 
formalism, one starts from a different lagrangian defined in a different function space, follows a different method and 
ends up in a different set of field equations. Bearing this in mind, one nevertheless sees the striking resemblance of 
the expressions obtained via the two variational procedures (Equations (p^|-p2[)). We shall now see this is not a mere 
formal coincidence but that the reason lies in the fact that all our test Lagrangians are diffeomorphism covariant. 

All previous cases can be considered as specializations derived from a very general Lagrangian n— form constructed 
locally as follows. 



L 



[gab, /^a.gab, ■ • • , (ai . . . /^a^ )5a6, l/', ^a^^^ ■ • • : ^(ai ... A^a,) V", , (23) 



that is, L is a function of dynamical fields g, ip and also of other fields collectively referred to as 7. Referring to g and 
ip as (f), L is called /—covariant, / € Diff(M), or simply diffeomorphism covariant if, 

Lif*m^f*m, (24) 

where /* denotes the induced action of / on the fields (p. Note that this definition excludes the action of /* on /7 or 
the other fields 7. It immediately follows that our test Lagrangians considered previously satisfy the above definition 
and as a result are diffeomorphism covariant. 

It is a very interesting result, first shown by Iyer and Wald (Rcf. |2l|) that if L in (^3|) is diffeomorphism covariant, 
then L can be re-expressed in the form 

L — L [gab, Rb (ai . . ai ) ^) , (25) 

where V denotes the Levi-Civita connection of gab, Rabcd denotes the Riemann curvature of gab and m = max(fc — 
2,^ — 2). Notice that everything is expressed in terms of the Levi-Civita connection of the metric tensor and also that 
all other fields 7 are absent. 

Applying the Iyer- Wald theorem in our test Lagrangians we immediately see that we could have re-expressed them 
from the beginning in a form that involves only the Levi-Civita connection and not the original arbitrary connection 
and vary them to obtain the corresponding 'Hilbert' equations. As we showed above, we arrived at this result 
by treating the associated Lagrangians as different (indeed they are!) according to whether or not they involved an 
arbitrary (symmetric) or a Levi-Civita connection. 

III. CONFORMAL STRUCTURE AND WEYL GEOMETRY 

For the more general nonlinear Lagrangians of the form f{q) where q = R, RabR"''', or Rabcd.R'^^'^'^ where / is an 
arbitrary smooth function considered in the previous section, the field equations obtained by the metric-connection 
formalism are of second order while the corresponding ones obtained via the usual metric variation are of fourth 
order. This result sounds very interesting since it could perhaps lead to an alternative way to 'cast' the field equations 
of these theories in a more tractable, reduced form than the one that is usually used for this purpose namely, the 
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conformal equivalence theorem (Re f. ^21). In this way, certain interpretational issues related to the question of the 
physicality of the two metrics (Ref. |23[p4|) associated with the conformal transformation would perhaps be avoided. 
As discussed in the previous section, the constraint (03) for Weyl geometry becomes 



ic (5,r,A) = A 



ab 



^Ib ~ ilb} ~ 2-^'''" {Qagmb + Qbgarn ~ Qjngab) 



(26) 



In order to examine the consequences of the Weyl constraint (^6|) we now apply the constrained first order formalism 
to the Lagrangian L = f (R) . Variation with respect to the Lagrange multipliers recovers the expression (|lj) of the 
Weyl connection. Variation with respect to the metric yields the equations 

f'R{ab) - If gab + Bab - 0, (27) 

where Bab is defined by (|l^). Variation with respect to the connection yields the explicit form of the Lagrange 
multipliers, namely, 

K"' = [Q^r - ^V') + b^c (Q'f - F'r) - g"' (Qcf ~ PJ') ■ (28) 



Substituting back this last result into equation ( p^ ) we find, 

Bab = 2Q(aPb)f' - P(aPb)f' + f P (aQb) - f'QaQb - gab (2Q™^™/' - f - Pf + f'P'^Qm) ■ (29) 



Inserting this result into equation (^7|) we obtain the full field equations for the Lagrangian L — f (i?) in the framework 
of Weyl geometry, namely, 

f'Riab) - If gab - PaPbf + gabPf = Mab, (30) 

where Mab is defined by 

Mab = ~2Q^aFb)f' - f'N(aQb) + f'QaQb + ffah {^Q.r.P"'!' - f + f'P^Qm) • (31) 

It is interesting to note that the degenerate case Qa = corresponds to the usual field equations obtained by the 
Hilbert variation in the framework of Riemann geometry, namely, 

f'Rab - If gab - VaVbf + gab^f = 0. 

It is known that these equations are conformally equivalent to Einstein equations with a self-interacting scalar field 
as the matter source (Ref. [2^ ). In what follows, we generalize this property of the f (R) field equations in Weyl 
geometry. To this end, we define the metric g conformally related to g with /' as the conformal factor. Under a 
conformal transformation, the Weyl vectorfield transforms as 

and the field equations (^0|) in the conformal frame read, 

f'R{ab} ~ Ijjgab - PaPbf +gabPf' = M^b, 

where ^ = /7, p = 1'^^ aP b = ifT^ P and M^b is given by 

Mab = f'QaQb - f'Pia Q b) - P aP bf + gab (^/'^"Q™ - f'Q^ + P/') • (32) 

Introducing the scalar field (p — In/' and the potential V{(p) in the 'usual' form (Ref. ^ 

vi^) = ^{f{R)r'iRf'iR)~f{R)), (33) 
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we find that the field equations take the final form, 

Gab = M?b - 9abV (v?) , (34) 

where, 

Gab = R{ab) — 2^9ab 

and 

M^b = QaQb ~ '^(aQb) + dab (^"Q™ " Q') • 

The field equations (^) could also be obtained from the corresponding conformally transformed Lagrangian using the 
constrained formalism. In that case, the equation of motion of the scalar field ip derived upon variation with respect 
to if implies that the potential V {ip) is constant. 

Formally, the field equations (p4) look like Einstein equations with a cosmological constant and a source term 
depending on the field Qa- However, they cannot be identified as such unless the geometry is Riemannian, 
i.e. Qa = 0. This will be the case only if the original Weyl vectorfield is a gradient, Qa — /^a^- Then it can be 
gauged away by the conformal transformation Ijab — (exp $) gab and therefore the original space is not a general Weyl 
space but a Riemann space with an undetermined gauge (Ref. p5| ). This was the case of the unconstrained method 
applied to the Lagrangian L = f (i?), wherein the Weyl vector was deduced using equation (^) and turned out to be 
Qa = 'Pa (111/')- (I'^ I^^f. |2^, this pecuharity is used in order to find out a subclass of theories based on a general 
D-dimensional dilaton gravity action, for which both unconstrained method and Hilbert variation yield dynamically 
equivalent systems.) This fact shows that unconstrained variations cannot deal with a general Weyl geometry and 
correspond to a degenerate case of the constrained method — the field equations obtained from the former can be 
recovered only by choosing specific forms of the Weyl vectorfield (Ref. p6|). 



IV. DISCUSSION 

The results obtained in Section ^ and Section ^ have the interpretation that a consistent way to investigate gen- 
eralized theories of gravity without imposing from the beginning that the geometry is Riemannian, is the constrained 
first order formalism. Applications to quadratic and / (i?) Lagrangians in the framework of Riemannian and Weyl 
geometry reveal that unconstrained variational methods are degenerate cases corresponding to a particular gauge and 
that the usual conformal structure can be recovered in the limit of vanishing Weyl vector. 

The generalization of the result stated above to include arbitrary connections with torsion can be an interesting 
exercise. The physical interpretation of the source term (equations (|30| ) and (|3|)) is closely related to the choice of the 
Weyl vectorfield Q. However, it cannot be interpreted as a genuine stress-energy tensor in general since, for instance, 
choosing Q to be a unit timelike, hypersurface-orthogonal vectorfield, the sign of MabQ'^Q^ depends on the signs of 
/' (i?) and the 'expansion' p aQ"" ■ 



The generalization of the conformal equivalence theorem presented in Section III opens the way to analysing 
cosmology in the framework of these Weyl / (i?) theories by methods such as those used in the traditional Riemannian 
case. The first steps in such a program may be as follows (Ref. |27| ): 

(A) Analyse the structure and properties of Friedmann cosmologies, find their singularity structure and examine 
the possibility of inflation. 

(B) Consider the past and future asymptotic states of Bianchi cosmologies. Examine isotropization and recoUapse 
conjectures in such universes. Look for chaotic behaviour in the Bianchi VIII and IX spacetimes. 

(C) Formulate and prove singularity theorems in this framework. This will differ from the analysis in the Riemannian 
case (cf. Ref. 22) because of the presence of the source term Mah- 



All the problems discussed above can be tackled by leaving the conformal Weyl vectorfield Qa undetermined while 
setting it to zero at the end will lead to detailed comparisons with the results already known in the Riemannian case. 
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